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A KINETIC  NETWORK  MODEL  FOR  NONLINEAR  VISCOELASTIC  FLOW 
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PROPERTIES  OF  ENTANGLED  MONODISPERSE  POLYMERS 

High  polymers  in  concentrated  solutions  or  in  the  melt 
exhibit  shear-rate-dependent  viscosity  and  normal  stress  dif- 
ferences. In  transient  flows  upon  inception  of  steady  shear,  the 
stress  growth  function  of  the  polymer  overshoots  through  a maxi- 
mum emd  then  decays  to  the  steady-state  value  asyn^totically 
for  large  shear  rates.  However,  for  small  shear  rates,  it  only 
increases  monotonically  to  its  limiting  steady-state  value.  Upon 
cessation  of  steady  shear,  the  stress  relaxation  function  de- 
creases monotonically  to  zero,  the  rate  of  which  is  higher  for 
larger  shear  rates.  These  rheological  properties  have  been  in- 
terpreted by  the  shear-induced  chemges  in  the  steady-state  entem- 
glement  density  in.  the  polymer  (1) . 

Lodge  explained  the  nonlinear  viscoelastic  flow  properties 
by  a "rubberlike  liquid"  model,  which  however  still  predicted  a 
shear-rate-independent  viscosity  and  primary  normal  stress  co- 
efficient. Many  modified  theories  (3-6)  have  since  been  pro- 
posed to  refine  this  model.  Most  of  them  en^loy  en^irical  ex- 
pressions for  the  memory  functions  to  give  the  rate-dependent 
properties. In  this  work,  we  shall  propose  a kinetic  network  model 
which  appears  to  predict  the  observed  nonlinear  viscoelastic  flow 
properties  of  entangled  polymer  systems. 

The  basic  molecular  roech2mism  of  this  model  attributes  the 
rate -dependence  of  viscosity  and  normal  stress  difference  to  the 
decrease  in  entemglement  density  with  increasing  shear  rate.  In 
a flow  field,  entanglements  are  being  formed  and  disengaged. 
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constantly.  The  entanglement  loss  process  is  caused  primarily 
by  the  imposed  shear  and  is  assumed  to  be  proportional  to  the 
shear  rate  to  a power  aj  ^ 


(1) 


where  n.  is  the  entanglement  loss/rate,  k.  is  the  rate  constant, 


, ^ , 

Y is  the  shear  rate,  n is  the /fcurrent  number  of  entemglements  along 


a representative  chain,  and^a(£l)  accounts  for  the  elastic  nature 
of  the  polymer.  The  entanglement  creation  process,  driven  by 


thermal  diffusion,  is  assumed  to  be  independent  of  shear  rate.. 


In  other  words. 


He  = ( ^ ^"o" 


(2) 


where  n^  is  the  entauiglement  creation  rate,  is  the  rate 


constant,  X is  a characteristic  time  determined  by  the  rate  of 
diffusion,  and  is  thus  a function  of  molecular  weight  and 


temperature,  a is  the  elasticity  parauneter,  emd  n^  is  the 


saturation  number  of  entanglements  along  a representative  chain. 


(n^-  n)  is  the  number  of  vacant  sites  where  entanglements  have 


been  lost  amd  can  be  likely  regenerated. 

At  steady  state,  the  creation  rate  equals  the  loss  rate 


(n 


n.)  auid  the  entanglement  density  stays  constant.  Equating 


c 

Eqs.  1 amd  2,  we  have 


P = n/n  = 


1 + 


(kj^Ac)  (Xy)' 


(3) 


Eq.  3 describes  the  probability  for  entamglements  to  exist  at 
different  shear  rates.  The  viscous  drag  from  the  entanglements 
is  believed  to  be  the  major  source  of  contribution  to  polymer 
viscosity  (1).  Therefore,  to  a good  approximation,  the  reduced 


viscosity  as  a function  of  shear  rate  is  also  given  by 


I 

i 


nil)  . 1 (4) 

"o  l+(kjA^)  (X4)“ 


where  is  the  zero- shear-rate  viscosity  of  the  polymeric  fluid. 

The  normal  stresses  come  from  the  elastic  energy  stored  in 
the  system,  emd  are  primarily  sustained  by  pairs  of  entanglements, 
because  the  elastic  deformation  energy  can  be  trapped  in  between 
these  slowly  moving  points.  The  primary  normal  stress  difference 
should  thus  be  proportional  to  the  square  of  the  probability  for 
individual  entanglements  to  remain  under  the  flow  situation,  which 
in  turn  corresponds  to  the  probability  for  pairs  of  entanglements 


to  exist  simultameously . In  other  words. 


^11  ~ ^ 
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where  9 is  a proportionality  constant  which  determines  the 
value  of  primary  normal  stress  difference.  The  reduced 
viscosity  and  primary  noinnal  stress  difference  as  fvinctions  of 
shear  rate  calculated  by  Eqs.  4 and  5 for  a equals  to  0.85  is 
given  in  Fig.  1.  The  curves  show  qualitatively  the  expected 
trend. 

Under  tremsient  flow  conditions,  such  as  inception  amd 
cessation  of  steady  shear,  the  entanglement  probability  P becomes 
a function  of  time  which  also  has  some  parametric  dependence  on 
shear  rate.  It  can  still  be  calculated  by  the  following  dif- 
ferential equation  governing  the  relationship  between  the  chamge 
of  entanglements  to  the  two  competing  rate  processes,  i.e.. 


n * He  - (6) 

where  n is  the  rate  of  change  of  enternglements . Inserting  the 
expressions  for  entanglement  creation  and  loss  rates.  Eg.  6 cem 

be  integrated  to  give  for  start“Up  experiments. 

. a 

k Y t 

P(t;Yo)  - P»(Yo^  + (1  - p^(Yp))exp(--— . — ) (7) 

^00  0 

where  P„(Yq)  is  the  steady-state  entemglement  probability;  emd 
for  stress  relaxation  experiments. 

• • ^c 

P(tn^)  = 1 - (1  - P„(Yo))exp(-  -|  t)  (8) 

X 

where  P„(y)  is  again  the  steady-state  entemglement  probability, 
in  this  case,  before  the  flow  is  stopped. 

Polymer  network  possesses  certain  elastic  characteristics. 
When  it  is  suddenly  deformed  upon  the  inception  of  steady  shear, 
the  shear  stress  bom  by  an  average  entemglement  first  increases 
with  time  rapidly  and  then  gradually  levels  off  towards  the 
steady-state  value.  A Maxwellian  model  is  assumed  to 
represent  this  behavior,  i.e.. 


o (t;Y-) 

; = 1 - exp(-  t/x) 


a (y.) 

00  Q 


where  0(t;Y  ) and  o (y  ) are  the  current  average  stress  emd 
0 0 

the  final  steady-state  stress  carried  by  individual  entanglements 
respectively,  and  r is  the  system  time  constmt.  The  macro- 


scopic stress  exerted  by  the  fluid  is  given  by  the  product  of 
the  average  microstress  supported  by  a remaining  entanglement 


emd  the  number  of  such  stress-bearing  entanglements.  The 


normalized  shear  stress  growth  functions  thus  calculated  are  pre- 
sented in  Fig.  2.  The  curves  qualitatively  reproduce  the 
expected  behavior.  Quemtitative  comparison  of  this  theory  with 
experimental  data  is  now  being  investigated,  the  results  of 
which  will  be  reported  in  forthcoming  publications. 
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? CAPTIONS 

1.  Theoretical  curves  of  log  reduced  viscosity  and  primary 

^ normal  stress  difference  vs  log  reduced  shear  rate,  where  X 

is  the  characteristic  time,  and  a is  the  elasticity  parameter. 
(The  ratio  is  taken  as  1,  emd  a = 0.85) 

2.  Normalized  stress  growth  function  for  three  different 
shear  rates  calculated  by  the  theory.  The  parameters  used 

I are  a » 0.85,  = 0.3,  kj^/k^,!®  = 1,  t = 0.5. 
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